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IRAS 16293-2422: Evidence for Infall onto a Counter-Rotating 

Protostellar Accretion Disk 

Anthony J. Reniijan^’^ and J. M. Hollis^ 

ABSTRACT 

We report high spatial resolution VLA observations of the low-mass star¬ 
forming region IRAS 16293-2422 using four molecular probes: ethyl cyanide 
(CH3CH2CN), methyl formate (CH3OCHO), formic acid (HCOOH), and the 
ground vibrational state of silicon monoxide (SiO). Ethyl cyanide emission has 
a spatial scale of ~20" and encompasses binary cores A and B as determined 
by continuum emission peaks. Surrounded by formic acid emission, methyl for¬ 
mate emission has a spatial scale of ~6"and is conhned to core B. SiO emission 
shows two velocity components with spatial scales less than 2" that map ~2" 
northeast of the A and B symmetry axis. The redshifted SiO is ~2" northwest 
of blueshifted SiO along a position angle of ~135° which is approximately par¬ 
allel to the A and B symmetry axis. We interpret the spatial position offset in 
red and blueshifted SiO emission as due to rotation of a protostellar accretion 
disk and we derive ~1.4 Mq interior to the SiO emission. In the same vicinity, 
Mundy et ah (1986) also concluded rotation of a nearly edge-on disk from OVRO 
observations of much stronger and ubiquitous ^^CO emission but the direction of 
rotation is opposite to the SiO emission Endings. Taken together, SiO and ^^CO 
data suggest evidence for a counter-rotating disk. Moreover, archival BIMA array 
^^CO data show an inverse P Cygni profile with the strongest absorption in close 
proximity to the SiO emission, indicating unambiguous material infall toward 
the counter-rotating protostellar disk at a new source location within the IRAS 
16293-2422 complex. The details of these observations and our interpretations 
are discussed. 

Subject headings: ISM: abundances - ISM: clouds - ISM: individual (IRAS 16293- 
2422) - ISM: molecules - radio lines: ISM 
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1. INTRODUCTION 

IRAS 16293-2422 is a low-mass star forming region located in the p Ophiuchus cloud 
complex at a heliocentric distance of 160 pc. It contains an undetermined number of proto- 
stellar objects, high velocity outflows in the E-W and NE-SW directions (Stark et ah 2004), 
and the region is dominated by two radio continuum peaks designated cores A and B. These 
two principal cores are separated by and are often referred to as the “binary” system 
with a mass of 0.49 Mq for core A and 0.61 Mq for core B (Looney, Mundy & Welch 2000). 
Moreover, the gas and dust toward core A appears to be at a higher temperature than core 
B and most of the molecular emission lines of high energy transitions have been detected 
exclusively toward core A. The measured dust temperature toward core B is ~40 K (Mundy 
et ah 1986) and the measured temperature toward core A from many different molecular 
species is between 80 and 200 K (Chandler et ah 2005). 

Spectral observations with the IRAM 30-m radio telescope toward the low-mass pro- 
tostellar system IRAS 16293-2422 (Cazaux et ah 2003) demonstrated an extremely rich or¬ 
ganic inventory with abundant amounts of complex oxygen and nitrogen bearing molecules 
including formic acid (HCOOH), methyl formate (CII 3 OCHO), acetic acid (CH 3 COOH) and 
methyl cyanide (CH 3 CN) which are archetypal species often found in massive hot molecular 
cores (HMCs). Subsequent spatial observations by Kuan et al. (2004) with the Submillimeter 
Array (SMA) mapped emission from high energy transitions of CH 3 OCIIO, vibrationally ex¬ 
cited vinyl cyanide (CH 2 CHCN) and other large molecular species. Additionally, Bottinelli 
et al. (2004) mapped slightly lower energy transitions of CH 3 CN and CH 3 OCHO with the 
Plateau de Bure interferometer, noting that all the molecular emission is contained toward 
cores A and B and that no large molecule emission is seen associated with the molecular 
outflows. Recently, Chandler et al. (2005) detected large abundances of a number of dif¬ 
ferent organic molecular species primarily located toward core A from SMA observations. 
Moreover, from observations of sulfur monoxide (SO) absorbing against the background dust 
continuum near binary core B, Chandler et al. (2005) suggest that there must be material 
from an extended molecular envelope falling onto an embedded dust disk and that there may 
be a low-velocity outflow coming from an embedded protostar. 

The discovery of large oxygen and nitrogen bearing molecules in a low-mass star forming 
region is signihcant because the molecular complexity appears to be similar to high-mass 
star forming regions. Such complex molecules undoubtedly are incorporated into subsequent 
formation of protoplanetary bodies and would provide a rich organic inventory that could 
be supplied to early planetary systems similar to that expected to have occurred in our own 
pre-solar nebula. With this in mind, we conducted high resolution observations of IRAS 
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16293-2422 using the NRAO^ Very Large Array (VLA) using four molecular probes: ethyl 
cyanide (CH 3 CH 2 CN), methyl formate (CH 3 OCHO), formic acid (HCOOH) and the ground 
vibrational state of silicon monoxide (SiO) to determine the implication of their relative 
spatial locations. Surprisingly, the location of SiO is not co-spatial with either core A or B. 
This prompted us to process archival BIMA array data for an explanation. 


2. OBSERVATIONS 

Observations of IRAS 16293-2422 were conducted with the VLA in 2005 June 27 and 
28 in its BnC conhguration. Formic acid (HCOOH) and methyl formate (CH 3 OCHO) were 
observed simultaneously in one correlator setting and ethyl cyanide (CH 3 CH 2 CN) and silicon 
monoxide (SiO) in another correlator setting. Table 1 lists the molecule (col. [1]), transition 
quantum numbers (col. [2]), rest frequency (col. [3]), upper energy level E„ (col. [4]), line 
strength S (col. [5]), dipole moment Ha (col. [ 6 ]), and molecular parameter reference (col. 
[7]). The phase center for all observations was a= 16^32™'22^.8, 6 = -24"28'33" (J2000.0). 
The correlator was operated in the four intermediate-frequency (IF) normal mode, and the 
selection of a 6.25 MHz bandpass yielded 32 spectral channels in two pairs of oppositely 
polarized IFs giving a channel width of 195.313 kHz ( 1.35 km s ^ at 44 GHz). 

The compact radio source 1625-254 was used as the primary phase calibrator. The 
absolute flux density scale was determined from observations of 1331-1-305 (3c286) whose 
average flux density from the 2 observing bands was 1.44 Jy as given by the SETJY routine 
in AIPS. The averaged bootstrapped flux density of 1625-254 was 0.6 Jy. Antenna gains were 
derived from 1625-254 observations at 5 minute intervals. After ~1 hr, a pointing calibration 
was performed to ensure pointing accuracy. All data were combined and calibrated using the 
AIPS data reduction package of NRAO. The calibrated u — v datasets were then imported 
into MIRIAD (Sault et al. 1995) for imaging and deconvolution. Figures 1-4 show the images 
of the emission from CH 3 CH 2 CN, CH 3 OCHO, HCOOH, and SiO, respectively toward IRAS 
16293-2422. 

The naturally weighted synthesized beamwidth resulting from the BnC conhguration 
at 44 GHz was 0."4x0."3 and the largest angular size that can be measured in the BnC 
array is 43". In order to determine the location and morphology of the molecular emission 
toward 16293-2422, the calibrated u — v data were INVERTed using synthesized beams 
ranging from the naturally weighted synthesized beamwidth to a beamwidth of 20 "x 20 ". 


^The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated 
under cooperative agreement by Associated Universities, Inc. 
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This process indicated whether the bulk of the emission from a given species was compact 
or extended and was utilized since: (1) emission from HCOOH, SiO (J = 1 — 0, n=0) and 
CH 3 CH 2 CN had not been previously mapped toward this region, and (2) no emission feature 
of low energy transitions (E,j<7.5 K) of large molecules had been imaged toward this region. 


3. RESULTS AND DISCUSSION 
3.1. CH 3 CH 2 CN 

Figure 1 shows a contour map of the 5 i 5 - 4 i 4 transition of CH 3 CH 2 CN superimposed 
on the 0.7 cm continuum emission (grayscale). The mapped CH 3 CH 2 CN emission region 
shows two blended emission peaks. Compared to the small naturally weighted beamwidth, 
a relatively large synthesized beamwidth of 10."0x7."5 is necessary to detect CH 3 CH 2 CN 
which is extended with respect to cores A and B and the oxygen bearing species (see § 3.2). 
Figure 1 also shows a banning smoothed emission spectrum, which was extracted from the 
entire region of mapped emission, with an rms noise level of ~7 mJy beam“^ shown on 
the left. The spectrum shows a relatively strong emission component at an LSR velocity of 
+7.7 km s“^ and a possible weaker component at +2.0 km s“^. The systemic LSR velocity 
of +3.9 km s“^ (Bottinelli et ah 2004) is shown as a dashed line. We determine the total 
beam-averaged CH 3 CH 2 CN column density, Nt, using the formula: 


Nt = 2.04 


rot)^ 




X 10^+m- 


( 1 ) 


which is described in Remijan et al. (2004) where Eu is the upper state energy level (K); T^ot 
is the rotational temperature (K); 6a and db are the major and minor axes of the synthesized 
beam (arcsec); J Aldv is the observed integrated line intensity (Jy beam“^ km s“^); i/ is the 
rest frequency (GHz); SijfA is the product of the transition line strength and the square of the 
electric dipole moment (Debye^) and Q{Trot) is the rotational partition function. A'r(EtCN) 
~ 5.2x10^^ cm“^, which assumes a rotational temperature of 54 K (Cazaux et al. 2003), is 
listed in Table 2. For a molecular hydrogen column density of 7.5x10^^ cm“^ (Cazaux et 
al. 2003), the CH 3 CH 2 CN fractional abundance, A(FtCN), is 6.9x10“®. These values are 
consistent to the column density of CH 3 CH 2 CN previously measured toward this source and 
is typical of the fractional abundance seen toward the Orion Compact Ridge (Cazaux et al. 
2003). 

The Figure 1 detection of extended emission with a spatial scale of ~20” in a low-energy 
transition of CH 3 CH 2 CN demonstrates that there is a large envelope of cold molecular gas 
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surrounding the compact binary cores A and B. Core A is reputed to be hot (Chandler 
et ah 2005) while core B is thought to be much cooler (Mundy et al. 1986). A warm gas 
phase chemistry or grain surface chemistry can lead to enhanced emission from nitrogen 
bearing molecules in hot molecular cores (Remijan et al. 2004; Mehringer & Snyder 1996). 
In previous observations of IRAS 16293-2422 which sample high energy level transitions at 
1 and 3 mm, enhanced emission of nitrogen bearing molecules have been detected toward 
cores A and B (e.g., CH 3 CN by Bottinelli et al. 2004; HCN, HC 3 N, and CH 2 CHCN by Kuan 
et al. 2004) on a typical spatial scale of ~5". Thus, our low-energy CH 3 CH 2 CN observations 
are complementary to these high-energy observations of other nitrogen bearing molecules, 
suggesting both increasing temperature and density gradients toward cores A and B. 


3.2. CH 3 OCHO and HCOOH 

High energy transition of a number of different molecules have emission centroids near 
core A (e.g., see Figure 9 of Chandler et al. 2005); in particular, CH 3 OCHO has high energy 
transitions that are associated with core A (Chandler et al. 2005; Bottinelli et al. 2004; 
Kuan et al. 2004). The only large oxygen bearing molecule detected by Chandler et al. 
(2005) toward core B is a low energy transition of dimethyl ether (CH 3 OCH 3 ). Figures 2 
and 3 show maps of low-energy transitions of CH 3 OCHO and HCOOH, respectively, that are 
in proximity to core B. This is further evidence that core B is a much cooler protostellar core 
than core A. In Figure 2, the averaged emission from both the A and E states of CH 3 OCHO 
are mapped using a synthesized beamwidth of ~2."6xl."8. Figure 2 also shows an emission 
complex centered on the A and E states if one assumes an LSR velocity of ~-|-3.9 km s“^. This 
complex appears to be suffering from the competing effects of emission and self-absorption. 
In order to observe both states in the same bandpass, the average frequency (45.3966 GHz) 
of both states was used in the observations and this average frequency corresponds to an 
LSR velocity of -1-3.9 km s“^ in Eigure 2 while the A and E states appear at -1.6 and -1-9.4 
km s“^, respectively. Figure 3 shows a map of the low-energy transition of HCOOH that 
surrounds the peak of core B and also the CH 3 OCHO emission in Figure 2; no HCOOH 
emission was detected in the vicinity of core A. In Figure 3, the emission from HCOOH was 
mapped using a synthesized beamwidth of ~ 2 ." 6 xl." 8 ; similar to the beamwidth used to 
map the CH 4 OCHO emission. The LSR velocity of HCOOH is ~-|-3.9 km s“^ as indicated 
by the vertical dashed line in the spectrum also shown in Figure 3. 

The spatial distributions of CH 3 OCHO and HCOOH toward core B of IRAS 16293-2422 
compare favorably to those same spatial distributions observed toward the Orion Molecular 
Cloud (OMC-1) compact ridge (Hollis et al. 2003; Liu et al. 2002). For example, in the OMC- 
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1 compact ridge, Hollis et al. (2003) observed HCOOH emission surrounding CH 3 OCHO 
emission, and suggested that the HCOOH emission delineates the leading edge of a shock 
front as the outflow from source ”1” interacts with the quiescent ambient gas; the CH 3 OCHO 
emission, which is closer to source ”1”, then represents the post shock gas. If this scenario is 
correct for IRAS 16293-2422, then the corresponding outflow source associated with Figures 

2 and 3 is likely core B itself. Indeed, Chandler et al. (2005) report evidence for material 
infall toward and also low-velocity outflow from core B; both conditions undoubtedly result in 
shock phenomena. Toward core A, enhanced abundances of large molecules provide indirect 
evidence for shocks to explain the enhancements (Chandler et al. 2005). The molecular 
differences between cores A and B are most likely explained by temperature differences and 
different protostellar ages. In any case, the presence of shocks seemingly leads to abundance 
enhancements in oxygen-bearing molecules, including large interstellar aldehydes located in 
high-mass star forming regions like Sgr B2 (e.g., see Hollis et al. 2004a, b). 

Because the CH 3 OCHO emission appears to be self-absorbed, we can only give a lower 
limit to the column density over the entire emission complex. Since we have observed only low 
energy transitions in this work, a conservative rotational temperature of 40 K (Mundy et al. 
1986) is assumed which yields ^^’(MeF) >1.9x10^® cm“^. Assuming an H 2 column density 
of 1.6x10^^ cm“^ toward core B (Kuan et al. 2004), the CH 3 OCHO fractional abundance is 
X(MeF) >1.2x10“®. For the same rotational temperature, AA(HCOOH) ~ 4.1x10^^ cm“^ 
and the HCOOH fractional abundance is A(HCOOH) ~ 2.5x10“® (see Table 2). 


3.3. SiO 

Figure 4 shows redshifted and blueshifted SiO (J = 1 — 0, v=0) emission contours 
mapped separately over the 0.7 cm continuum (grayscale). These two components of the SiO 
emission are seen at LSR velocities of -1-2.3 km s“^ and -1-7.9 km s“^ which are approximately 
symmetrical around the systemic LSR velocity of 3.9 km s“^. The peak intensity maps, which 
are reminiscent of SiO maser spots, are shown at these velocities. The resulting synthesized 
beamwidth of ~l."5x0."8 is shown at the bottom left of the map. Also shown in Figure 4 
is a spectrum which is averaged over the SiO emitting regions and banning smoothed over 3 
channels; the rms noise level is ~3 mJy beam“^ shown to the left. 

This is the hrst detection of the J = 1 — 0, n=0 transition of SiO toward IRAS 16293- 
2422 and the hrst map of the compact emission from SiO in the vicinity of cores A and B. 
Previous surveys by Blake et al. (1994) and Ceccarelli et al. (2000) detected weak emission 
features of higher energy transitions of n=0 SiO using single element radio telescopes. From 
a chemical model of the IRAS 16293-2422 region, Ceccarelli et al. (2000) concluded that 
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an enhanced abundance of SiO was necessary to account for the excess emission from the 
higher J level transitions and that the emission must be coming from a compact (<3") region. 
Finally, their model calculated an infall velocity of SiO of 2.8 km s“^ which was consistent 
with their observed ~5 km s“^ linewidths. 

From our high resolution VLA observations, we hnd that the blue and redshifted SiO 
emission peaks have a spatial separation of and the center of the SiO emission centroid is 
located at a= 16^32^227875, 5 = -24'’28'32".48 (J2000.0). This location is within 0."6 of the 
pointing position and location of the ^^CO emission centroid taken with the Owens Valley 
Radio Observatory (OVRO) millimeter array (Mundy et ah 1986) and the peak position of 
CS emission taken with the IRAM 30-m (Menten et ah 1987). From the red and blueshifted 
emission peaks, Mundy et al. (1986) inferred a rotational velocity of 2 km s“^ and a 2 Mq 
rotating disk of material internal to the ^^CO emission. From the CS data, Menten et al. 
(1987) inferred a rotational velocity of 1.1 km s“^ from lower resolution observations a 4 
Mq rotating disk of material. From the high resolution SiO data, we determine a rotational 
velocity of ~2.8 km s“^ and for a IRAS 16293-2422 distance of 160 pc, the mass interior 
to the SiO emission is ~1.4M0. However, assuming the SiO emission is associated with 
rotation, the compact disk containing the SiO emission is counter-rotating with respect to 
the previously reported ^^CO and CS emission. From archival BIMA array data taken of 
the J = 2 — 1 transition of ^^CO at 230.538 GHz, we hnd both emission and absorption 
toward IRAS 16293-2422. Moreover, the spectrum of the ^^CO data clearly show an inverse 
P Cygni prohle with the strongest absorption in close proximity to the SiO emission (Figure 
5), indicating unambiguously, material infalling toward the molecular disk. Chandler et al. 
(2005) also suggest there is evidence for infall toward core B from the absorption prohle of 
SO taken with the SMA. 

It is important to note that we have only observed small-scale, blue-shifted emission and 
small-scale, red-shifted absorption of ^^CO in the vicinity of cores A and B (see the image in 
Figure 5); these data possess the characteristic signature of infalling matter (i.e., an inverse P 
Cygni prohle shown in Figure 5). While the ^^CO emission encompasses cores A and B, the 
more compact absorption occurs at a discrete position ~3."5 north of core A and ~2."5 west 
of core B (see Figure 4). On the other hand, large-scale ^^CO has been observed in red-shifted 
emission and blue-shifted emission by other investigators who reasonably interpret these 
results as outhow whose origin has traditionally been assigned to core A. Thus, ^^CO on two 
diherent spatial scales indicate two diherent phenomena: one of infall on a small-scale toward 
a newly detected protostellar accretion disk whose position is precisely determined and one 
of large-scale outhow whose exact origin may be difficult to determine. The location of the 
accretion disk that we have detected seems to account for our small-scale SiO observations 
and the somewhat larger-scale ^^CO observations of Mundy et al. (1986) and CS observations 



of Menten et al. (1987) only if the nearly edge-on disk is counter-rotating since the direction 
of rotation is opposite for SiO when compared to either or CS. Additional observational 
support for this scenario can be seen in the BIMA and OVRO array images of Schoier et 
ah (2005) who observed a number of low energy transitions (i.e., the J = 1 - 0 of HCO'*', 
HNC, and N 2 H'^ and the J = 2 — 1 of and SiO). In these images, for example, the 
J = 2 — 1 transition of SiO suggests rotation in the opposite direction to the J = 1 — 0 
transition of HCO'*'. Compared to the VLA data presented in this work, these BIMA and 
OVRO data have an order of magnitude coarser spatial resolution that precludes an exact 
determination of the SiO emission locations that are apparent in the VLA map shown in 
Figure 4. Furthermore, SiO is a refractory molecule that can withstand high temperatures, 
and would likely exist closer to a forming protostellar source than most other molecules. 
Thus, the high resolution VLA observations have shown for the hrst time, evidence for a 
new protostellar core with a compact inner portion of the disk containing the SiO emission 
counter-rotating with respect to the outer portion of the disk containing ^^CO and CS. 
Counter-rotating cores have been reported in disks of galaxies before (e.g., Thomas et al. 
2005 and references therein), but this is the first indication of such behavior in a protostellar 
accretion disk. 


4. CONCLUSIONS 

To summarize, we hnd a low temperature envelope of molecular CH 3 CH 2 CN emission 
that is encompassing the IRAS 16293-2422 cores A and B. Furthermore, low frequency, low 
energy transitions of large oxygen bearing molecules (i.e. HCOOH and CH 3 OCHO) are seen 
toward core B where an outflow may be responsible for a low velocity shock that is leading 
to enhanced abundances of large oxygen bearing species. Also, the SiO emission shows two 
velocity components with spatial scales less than 2" and a velocity separation of ~5.6 km 
s“^. We interpret the spatial position offset in red and blueshifted SiO emission as due to the 
rotation of a protostellar accretion disk and we derive ~1.4 Mq interior to the SiO emission. 

While it could be that the SiO complex we observe is due to outflow, evidence is ac¬ 
cumulating to the contrary that rotation is the explanation. Such evidence is best seen 
on a small scale. The SiO observations in the present work show two velocity components 
each with linewidths of ~5 km s“^, suggesting rotation of a disk seen edge-on since the two 
components are not cospatial, are symmetric about an infall location, and have a velocity 
gradient perpendicular to large-scale CO outflow. Moreover, images of other molecules on a 
small-scale also show that their velocity gradients are orthogonal to the direction of outflow 
seen on a large scale (e.g., Schoier et al. 2005). In this work, we have been very careful to 



match the synthesized beam to the source emission or absorption to pin-point the location 
of a new source of infall ~3."5 north of core source A and ~2."5 west of core source B. We 
have done this by two independent sets of observations (i.e., the VLA in SiO emission and 
the BIMA array in absorption). While there are no previous reports of continuum 

emission at this location, an inverse P Cygni prohle is absolute proof of absorption against 
background continuum, and our BIMA data show that the location is smaller than the syn¬ 
thesized beam (i.e., 3.2"x0.8"). Thus, the continuum source is undoubtedly small (e.g. a 
disk seen edge-on) and may eventually be detected by an interferometer at high resolution 
with enough integration time. 

However, the compact disk containing the SiO emission is counter-rotating with respect 
to the and CS emission seen by other investigators. This is the hrst report of evidence 
for a counter-rotating accretion disk toward a low-mass protostellar complex. Moreover, 
archival BIMA array data show an inverse P Cygni prohle with the strongest absorption 
in close proximity to the SiO emission, indicating unambiguously, material infalling toward 
the counter-rotating protostellar disk. The infall location is compelling evidence for a new 
protostellar source within the IRAS 16293-2422 complex. 

We thank F. Schoier and C. Ceccarelli for valuable comments on this work, and R. 
Crutcher for permission to utilize the BIMA array data archive to explore ^^CO data toward 
IRAS 16293-2422. J.M.H. gratefully acknowledges research support from H.A. Thronson 
while assigned to the NASA Science Mission Directorate. 
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FIGURE CAPTIONS 

Fig. 1.— Spectrum and map of the 5 i 5 - 4 i 4 transition of GH 3 GH 2 GN. The spectrum, which 
was taken from the entire region of mapped emission, is banning smoothed over 3 channels. 
The rms noise level of ~7 mJy beam“^ is shown at the left. The systemic LSR velocity of 
+3.9 km s“^ is shown as a vertical dashed line. The contour map of the 5 i 5 - 4 i 4 transition of 
GH 3 GH 2 GN is shown superimposed on 0.7 cm continuum emission (grayscale). The contour 
levels are -0.045, 0.045, 0.060, 0.075 and 0.09 Jy beam“^ (0.045 Jy beam“^ = 5 cr). The 
synthesized beamwidth is shown at the bottom left of the map. 

Fig. 2.— Spectrum and map of the 5 i 5 - 4 i 4 A and E state transitions of GH 3 OGHO. The 
spectrum, which was taken from the entire region of mapped emission, is banning smoothed 
over 3 channels. The rms noise level of ~2 mJy beam“^ is shown at the left. The average 
frequency of the A and E state transitions appears at an LSR velocity of +3.9 km s“^. 
The contour map of the GH 3 OGHO emission is shown superimposed on 0.7 cm continuum 
emission (grayscale). The contour levels are -0.008, -0.010, 0.010, 0.012, 0.014 and 0.016 Jy 
beam“^ (0.010 Jy beam“^ = 5 cr). The synthesized beamwidth is shown at the bottom left 
of the map. 

Fig. 3.— Spectrum and map of the 2 o 2 -loi transition of HGOOH. The spectrum, which was 
taken from the southernmost region of mapped emission toward core B, is banning smoothed 
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over 3 channels. The rms noise level of mJy beam”^ is shown at the left. The spectrum 
shows the HCOOH emission is centered at the systemic LSR velocity of +3.9 km s“^ (vertical 
dashed line). The contour map of the HCOOH emission is shown superimposed on 0.7 cm 
continuum emission (grayscale). The contour levels are -0.007, -0.009, 0.011, 0.014, 0.016 
and 0.018 Jy beam“^ (0.011 Jy beam“^ = 5 a). The synthesized beamwidth is shown at the 
bottom left of the map. 

Fig. 4.— Spectrum and maps of the J=l-0, n=0 transition of SiO. The spectrum of SiO was 
averaged over the emission regions and is banning smoothed over 3 channels. The spectrum 
shows a redshifted and blueshifted components of SiO emission approximately centered at 
the systemic LSR velocity of +3.9 km s“^ (vertical dashed line). The 1 a rms noise level of 
~2 mJy beam“^ is shown on the left. The contour map of the SiO emission is superimposed 
on the 0.7 cm continuum emission (grayscale) and shows the location of the redshifted and 
blueshifted emission peaks. The contour levels are -0.008, -0.010, 0.010, 0.011, 0.012, 0.013 
and 0.014 Jy beam“^ (0.010 Jy beam“^ = 5 a). The synthesized beamwidth is shown at the 
bottom left of the map. Also shown in the map is the location of the absorption trough of 
the inverse P Cygni prohle (hlled circle) of the emission shown in Figure 5. 

Fig. 5.— Map and spectrum of the J=2-l transition of at 230.538 GHz taken with 
the BIMA array. The contour map of the emission is superimposed on the 1 mm 

continuum emission (grayscale) and shows the location of the redshifted absorption and 
blueshifted emission. The contour levels are -5.50, -5.00, -4.50, -4.00, -3.50, 2.25, 3.00, 3.75, 
4.50, 5.25 and 6.00 Jy beam“^. The synthesized beamwidth of 3."2x0."8 is shown at the 
bottom left of the map. The spectrum of the J=2-l transition of at 230.538 GHz was 
taken at a location surrounding the absorption and emission peaks and has a 1 cr rms noise 
level of ~1 Jy beam“^ shown on the left of the spectrum. The systemic LSR velocity of +3.9 
km s“^ is shown as a vertical dashed line. 



Table 1. Molecular Line Parameters for VLA Observations 



Transition 

Frequency®" 

Em 


l^a 


Molecule 

+ 

1 

'+ 

1 

(MHz) 

(K) 

S 

(D) 

Reference 

(1) 

(2) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

SiO 

1 

0 

0 

43,423.864(22) 

2.085 

1.00 

3.10 

b 

CH 3 CH 2 CN 

5i5-4i4 

43,516.200(8) 

7.380 

4.80 

3.85 

c 

HCOOH 

2o2-1o1 

44,911.750(50) 

3.238 

2.00 

1.39 

d 

CH 3 OCHO 

4i4-3i3 E 

45,395.755(16) 

6.187 

3.80 

1.63 

e 


4i4-3i3 a 

45,397.394(16) 

6.187 

3.80 

1.63 

e 


^Uncertainties in parentheses refer to the least signihcant digit and are 2 a 
values. 

'^(Lovas & Krupenie 1974) 

"(Lovas 1982) 

'^Willemot et ah 1980 
®Oesterling et ah 1999 


Table 2. Molecular Column Densities 



vlsr Range 

/Idv 

T 

rot 

Nt 

X 

Molecule 

(km s“^) 

(Jy bm“^ km s“^) 

(K) 

(xlO^^ cm“^) 

0 

1 

CD 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

CH3CH2CN 

+0.0 to +10.0 

0.34(9) 

54 

5.2(1.4) 

6.9(1.4) 

CH3OCHO 

-5.0 to +16.0 

>0.22 

40 

>187 

>11.7 

HCOOH 

+0.0 to +8.0 

0.07(3) 

40 

41(17) 

2.5(1.7) 


Note. — SiO is likely masering and is therefore not included in this table (see text). 
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